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The magnetization curves as a response of sweeping mag- 
netic field in the thermal environment are investigated using 
the quantum master equation. In a slow velocity region where 
the system almost behaves adiabatically, the magnetic plateau 
appears which has been observed in the recent experiment of 
Vi5 [ Phys. Rev. Lett. (2000)]. We investigate its mechanism 
and propose that this phenomenon is quite universal in the 
quasi-adiabatic transition with small inflow of the heat, and 
we call it 'Magnetic Foehn Effect'. We observe the crossover 
between this mechanism and the Landau-Zener-Stiickelberg 
mechanism changing the velocity. Some experiment is pro- 
posed to clarify the inherent mechanism of this effect. 



PACS number: 75.40.Gb,76.20.-fq ^ 

Properties of quantum dynamics have been studied ex- 
tensively for nanoscale magnets and also in microscopic 
system with nanostructure . There the nonadiabatic 
transition plays important roles. Landau Zener 
and Stiikelberg (LZS) derived the well-known nonadi- 
abatic transition probability for a two-level system with 
a sweeping field. It depends on the sweeping velocity 
and energy gap. Although the LZS formula is derived 
in two-level systems, it can also describe nonadiabatic 
transitions at the avoided crossings in the uniaxial mag- 
nets 0,1 . Therefore it can be widely applied to analyze 
phenomena related to the nonadiabatic transition in var- 
ious materials However real experiments are al- 
ways done in a thermal environment. Hence the studies 
on the nonadiabatic transition with an environment are 
quite crucial for further understanding of time-dependent 
phenomena in magnetic systems There universal as- 
pects of thermal effect independent of details of reservoirs 
are very important. As an example of such universal as- 
pects, we have studied the magnetization process in uni- 
axial molecule magnets such as Mni2 and Fcg at very 
low temperatures in thermal environment pO[ |. There 
the step- wise magnetization process is observed, which is 
due to the nonadiabatic transition and a fast damping 
process to the ground state. We call it apparent (or de- 
ceptive) nonadiabatic process. We also found that the 
transition probability of purely quantum mechanics is 
deducible from this process. This mechanism does not 
depend on the detailed structure of reservoir apid fluctu- 
ation of noise at the resonant field 

In this Letter, we propose another universal qualitative 
aspect of the thermal effect on the quasi-adiabatic tran- 
sition where the LZS probability is almost one. We un- 
derstand that the effect is quite general in systems which 
behave almost adiabatically in the dissipative environ- 



ment. The present study is directly related to the recent 
experiment for the molecule V15 which is effectively re- 

This molecule has 15 



garded as two level system 1 1 1 
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atoms of V but they divided into subclusters of 6,3, and 6 
spins. The subcluster of 6 spins forms a singlet state and 
contributes little to the magnetization and only the clus- 
ter of 3 spins mainly contributes to the magnetization. 
Therefore the model can be regarded as a two level sys- 
tem. This molecule is a very simple system and we may 
expect to see the LZS process clearly Q-^. However it 
was observed that the scattered population (i.e., that at 
the excited level) decreases when the sweeping velocity 
becomes faster |lj]. This is opposite to what we expect 
in view point of the LZS mechanism where the probabil- 
ity of the adiabatic transition should decrease for faster 
change of the field and the population of the excited state 
should increase. 

Chiorescu et al. jl^] explained this behavior in the view 
of the phonon bottleneck effect which means a lack of 
phonon number which contributes to the excitation from 
the ground state at near resonant magnetic fields. In 
the experiment, the heat reserver has a double-structure, 
i.e., the spin system is attached to the phonon system 
of the crystal, and the phonon system is attached to the 
external reservoir which is the liquid He. The contact 
between the phonon system and the external reservoir is 
so week that in short time scale, thermal effect in the spin 
system is caused by only phonons. Due to this eff'ect, the 
population of the excited state does not increase enough, 
and saturates at some value, which causes a magnetic 
plateau. 

In this Letter, we show that plateau in the magneti- 
zation curve is also observed in the case that the spin 
system is connected with a single heat reserver with slow 
relaxation rate, and propose that this qualitative prop- 
erty is universal with regardless to the detailed structure 
of environment. We investigate a magnetization pro- 
cess for a sweeping field by making use of the quantum 
master equation which we have used in our studies of 
quantum dynamics in dissipative environments pO| , p^ . 
Thereby, we investigate the magnetization plateau for 
various sweeping velocities and temperatures. 

The Hamiltonian we shall consider is given by 



-H = H{t) 



(1) 



where H{t) is the sweeping field, H{t) = vt and F is 
the transverse field. Transverse field represents a term 
causing quantum fluctuation and does not commute with 
the magnetization Sz- This simple system is realized in 
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many cases, e.g., the isotropic anti- ferromagnetic Heisen- 
berg chain with odd number of spins has the doublet in 
the ground state. ActuaUy V15 is in this situation [T^ . 

For this system (|l|), the LZS transition probabihty is 
given by, 



Phzs{v) = 1 - cxp 



V 



(2) 



in the case of 5 = 1/2 spin system |y-|6[. Thus the nor- 
mahzed magnetization at i = cxd is given by. 



1-2PlzsM- 



(3) 



We should note that this expression of A/out is also exact 
for any values for S although (j^) is derived for S = ^. 

We introduce a thermal environment taking the 
phonon system as the bath, TCb = J^Lj^^Lbui, where b^^ 
and are the annihilation and creation boson opera- 
tors of the frequency iv. We adopt the spectral density 
of the boson bath /(w) in the form /(w) = /qw" (w > 
0), (w < 0), with a = 2. In the experimental situation 
in the magnetic molecules such as S' = 10 in Mni2 and 
Fcg, the hyperfine interaction and the dipole interaction 
are not negligible at very low temperatures [ p^lp^ ]. In 
the case of V15 , the phonon gives a dominant contribu- 
tion 0. 

In the case of phonon bath, we can derive an equation 
of motion of the reduced density matrix p tracing out the 
degree of freedom of the bath in the following form (the 
quantum master equation fl^ ): 



[X,Rp{t)] + [X,Rp{t)Y\ , (4) 



where X is a system operator through which the system 
and the bath couple with the constant A. The first term 
of the right-hand side describes the pure quantum dy- 
namics of the system while the second term represents ef- 
fects of environments at a temperature T{— (3^^). There 
R is defined as follows: 

{k\R\m) = CiEk - E,n)ni3{Ek - E„,){k\X\m), 

C{uj) = I{uj)- and n^(cj) = (e'^"' - l)-\ 

where \k) and \m) represent the eigenstates of Ti. with 
the eigenenergies Ek and Em, respectively. 

We simulate the evolution given by Eq. (||) for various 
sweeping velocities in 1 /2 spin system of (|^) . Throughout 
this Letter, we set h to be unity. From now on, we set 
parameters as F = 0.5, T = 1.0, and A = 0.001. In 
Fig. 1(a), we present the magnetization curves for fast 
sweeping rates, v = 0.1, 0.2, and 0.4. Here we clearly find 
that the magnetic plateau decreases when v increases, 
which is consistent with (^. On the other hand, in the 
case of slow sweeping rates we also find the magnetic 
plateau as shown in Fig. 1(b) although in these sweeping 
rates the LZS transition probabilities (|^) are almost one. 
Here we should note that the magnetic plateau increases 



when when v increases, which is an opposite property to 
the fast sweeping case. This is the same phenomenon as 
the experiment p2[. 
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FIG. 1. Magnetization process for (a)« — 0.1,0.2, 
' = 0.4, and (b) v = 0.001, v = 0.006, and v = 0.01 
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In Fig. 2 we show the nonmonotonic dependence of the 
plateau height Afout on velocities. Trivially when v is 
very large, Afout shows linear dependence on 1/v from 
Eq. (^. In the slow sweeping rate region with Plzs ~ 1, 
A'/out goes down. 
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FIG. 2. The plateau height, Mout(f), as a function of the 
sweeping rate 

Here we discuss the population in the excited level 
P22(.H) for various sweeping velocities. Time depen- 
dences of p22{H) are shown in Fig. 3, where the popu- 
lation of p22{H) changes as a function of the sweeping 
velocity. Because the matrix element of the operator R 
is proportional to A and A", (A = ^/IP^+T^) , the ther- 
malization rate 7(A, A) is proportional to these values as 
7(A,A) cx AA". 

Let us investigate the relation between the thermaliza- 
tion rate 7(A, A) and the sweeping velocity v. If v is much 
larger, i.e., v 3> j{H,X), then no thermal relaxation oc- 
curs. In this case, abrupt change in the distribution p22 
due to the nonadiabatic transition takes place at i? = 0. 
This behavior is demonstrated in the case of w = 2 x 10~^, 
and 4 x 10~^ in Fig. 3. This sudden change causes the 
magnetic jump as shown in Fig. 1(a). Here it should be 
noted that we see the precession which means that the 
state is still highly coherent. In this region of the veloc- 
ity, when the velocity becomes small, the plateau goes 
up because Plzs monotonically increases until a simple 
adiabatic magnetization curve appears due to Plzs 1- 

In the further slow velocity, thermalization process be- 
gins to take place, that is, P22{H) tends to relax to its 
equilibrium value. 
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This increase of p22{H) causes a magnetic plateau as 
shown for w = 2 X 10"^, and 8 x 10~^ in Fig. 3. In 
order to realize a visible plateau, thermalization rate 
should be small comparing with v. Actually, because of 
7(if, A) cx AA^, thermalization is very slow at around the 
resonant point due to small A. We also studied the sys- 
tem with a = 0, where we found that the plateau sustain 
for large values of H{t) and the shape of the magnetiza- 
tion process seems very different. 
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FIG. 3. Field dependence of p22'- (a) the equilibrium value 
and (b) p22 in the sweeping field 



We associate the mechanism of magnetic plateau with 
the well-known Foehn phenomenon in the meteorology. 
The air with the vapor climbs up the mountain getting 
colder adiabatically and next the rain ensues. At this 
moment the vapor gives the heat to the air as the latent 
heat. Then, the air alone goes down the mountain and 
the temperature of the air increases higher than the orig- 
inal one due to the inflow of the latent heat at the moun- 
tain. In the present magnetic system, 'climbing up the 
mountain' corresponds to 'sweeping of magnetic field to 
H = 0', 'the latent heat' to 'the heat from the phonon', 
and 'the increase of temperature of the air' to 'the in- 
crease of temperature of the spin'. After the plateau 
the magnetization relaxes to the equilibrium value due 
to cooling by the thermal bath, which corresponds to the 
cooling of the air by the land after hot air reaches to 
the ground. Thus these similarities lead us to call this 
magnetic phenomenon 'Magnetic Foehn effect'. We con- 
firmed that the Magnetic Foehn effect is also observed in 
the Hamiltonian (|l|) with larger S. The essential mech- 
anism of this phenomenon is inflow of heat during adi- 
abatic process. Because the mechanism is quite simple, 
we can say that the 'Magnetic Foehn effect' is a univer- 
sal phenomenon in the magnetic systems which behave 
almost adiabatically. 

Because the LZS transition occurs only in the vicinity 
of H — and at other points only relaxations due to the 
dissipation term occur, the time evolution process may 
be divided into three regions:p(t) ^Ud{t : 0) -Uhzs -i^diO ■ 
—oo)p{—oo), where Ud(t : s) is the time evolution opera- 
tor due to dissipation part from a time s to t, and Ui,zs 
represents the LZS scattering matrix which corresponds 
to the pure quantum dynamics in the first term in (^) 
JlTf . For slow velocity case i; <C 1, we can consider that 
the system behaves adiabatically and thus we put the 
evolution Z/^lzs to be unity. Thus the transfer between 
the levels occurs only by the dissipation Ud{t : s). If we 
write down the second term of explicitly, we have the 
following equations for the diagonal term 

P22 = -^^i2An/3(A) C(A) [(e^^ + IV22 - 1] , (6) 

where X12 is the matrix elements of the operator X. In 
our simulation, it is expressed as X12 = I , ' for X = 
Sx + Sz ■ The term X12 depends on the choice of X, and 
is not universal. This equation leads the Magnetic 
Foehn effect. 

At the end of this Letter, let us discuss some experi- 
mental situations. As for the observation of the crossover 
from the Magnetic Foehn region to that of LZS region is 
difficult in V15, because the crossover from the Magnetic 
Foehn region to the LZS region locates at a very fast 
sweeping rate. For example, if F = O.IK, the crossover 
sweeping rate v is of order 10^ H/sec. Such fast change 
of the field is not easy to realize. However when the gap 
become small, the field at the crossover sweeping rate v 
becomes small. In Fig. 4, we show the dependence of the 
gap on the size and anisotropy for the system: 
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(7) 



where Si is an 5 = 1/2 spin and the open boundary con- 
dition is adopted. The gap becomes small with and 
decreases with N exponentially. Thus we expect that the 
crossover would be observed in some system and there we 
can obtain many informations of effects due to combina- 
tion of the quantum process and the thermal effect. 
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FIG. 4. Dependence of the gap on the size of the system A*' 
and the anisotropy A 

In Vi5, Chiorescu et al. attributes the Magnetic Foehn 
phenomenon to the lack of phonon as a mechanism of the 
insufficient supply of heat. In such a situation, we would 
propose an experiment to realize the adiabatic transi- 
tion of an isolated spin system. That is, we first sweep 
the field from the large negative value to the field (i/p) 
where the magnetic plateau (Mp) is observed. At this 
point the number of phonon is very small and supply 
of phonons from out side is expected to be slow. In 
this circumstance, if the field is swept in the opposite 



direction from H-p to 



-H-p, the spins at the lower spin 



can not be excited by the phonon because no phonon is 
available, and behaves pure quantum mechanically. In 
the experiment the LZS probability is almost one. The 
spins at the higher level may emit the phonon and re- 
lax to the lower level. But the emitted phonon will be 
used to excite the spin again because the population of 
the upper level is much smaller than that of equilibrium. 
Thus we expect that the magnetization simply changes 
the sign when iJp 



cess 



—Hp. In the iteration of this pro- 
Hp —Hp —f • • •, the magneti- 



zation would maintain the same amplitude for a while, 
-1 Mp -> -Mp ^ Mp ^ -Mp ■ ■ •, before the heat 
fiows in from the external environment and equilibrates 
the system. From this slow relaxation of magnetization 
we could know the relaxation rate between the phonon 
and the external bath. On the other hand, in the Mag- 
netic Foehn phenomenon due to slow relaxation but not 
short of phonon number, Mp relaxes with the thermal- 
ization rate. From Eq.(M), we can derive the relation of 



P22 for the iteration: — pi +P2P22 with ^22^ = 
and pi = {1 — Mp)/2. The sequence Mp"'' is given by 

|M(")| = l-2p^^^ where p^^^ = (l-Mp)(l-p^)/(l-p2). 
Here p2 is given by 



P2 = exp 



-A' 



dT A"(r) coth(/3A(T)/2) 



(8) 
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